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The model of a quantum-optical device for a conditional 
preparation of entangled states from input mixed states is 
presented. It is demonstrated that even thermal or pseudo- 
thermal radiation can be entangled in such a way, that Bell- 
inequalities are violated. 

PACS number(s): 03.65.Bz, 42.50.Dv 



I. INTRODUCTION 

In last decades, the phenomenon of entanglement be- 
tween two spatially separated photons was investigated 
both experimentally and theoretically mainly in order to 
show that quantum mechanics is not a local realistic the- 
ory 0| . As a counterpart to the particle-like behavior of 
photons, the entanglement of coherent states, that can 
be considered as the quantum analogue of deterministic 
light waves, was examined 0-0]. In both the cases, the 
entangled states, that were used to test Bell inequalities, 
were usually considered to be pure states. Recently, the 
entanglement of mixed states has been analysed to un- 
derstand, how the disorder influences on the amount of 
entanglement ||. 

In this report, we examine a new situation, when an 
entangling device prepares the entangled states of radia- 
tion from mixed states (thermal or pseudo-thermal light) 
at the input. Similarly to the idea presented in Ref. jjj, 
the entangling device can produce a four-mode entangled 
state with two mixed states and two vacuum states. It 
is shown, that even for very disordered states Bell in- 
equalities can strongly be violated. If there is a narrow 
frequency portion of thermal radiation in the input of 
the entangling device then Bell inequalities are violated 
when the frequency of radiation is "low" and the temper- 
ature of thermal source is "high" . For a pseudo-thermal 
radiation the violation of Bell inequalities is even more 
significant. In addition, the violation can be enhanced for 
both the cases of radiations, if a lot of different modes 
are entangled with vacuum state. Thus almost the max- 
imal Bell inequality violation can be achieved with such 
thermal states exhibiting a large entropy. 



II. PREPARATION OF MIXED ENTANGLED 
STATE 

We consider two separate systems A and B which con- 
sist locally of two modes Al, A2 and Bl and B2. All 
modes are initially uncntanglcd. We further assume that 
the density matrices of these four modes are diagonal in 
orthonormal Fock (number-state) bases {|n)} and that 
the modes A2 and B2 are in vacuum states, 



PA = y^Pn\n)Al(n\ ® 10)^42 (0|, 



PB = ^2r m \m) B i(m\ (g) |0) B2 (0|. 



(1) 



The density matrix of the total system has a factorized 
form /5i n = pa® Pb- Now, one can consider a condi- 
tional operation which enables to prepare the following 
entangled states (for n ^ or m ^ 0) 
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\lpnm) = (|^)ai|0)a2|0)bi|to)b2 

-|0)Ai|n) A 2|m) B1 |0) B2 ). 



(2) 



The entangling device prepares, for each m, n, the ana- 
logue of a singlet state, that was often employed to test 
Bell-type inequalities. Thus the initial density matrix p m 
is transformed into the form 



Pout — 

(1 - S n0 5 m0 )\tp nm ) (tp 

nrn \ • 



in n 



(3) 



where N = E nm Pn»ro(l - S n0 S m0 )] 1 = (1 - Po?*o) 1 - 

If there is at least one n > and one m > such 
that p n and r m ^ then state (||) is entangled. It 
can be proved in a very straightforward way using the so 
called transposition criterion . This criterion says that 
if operator p TB , obtained from p by partial transposition 
in subsystem B, is not positive the state p is entangled. 
Partial transposition of 

Pout = Pijklmnst \iAl3A2kBllB2){mAinA2SBltB2\ 

ijklrnnst 

in basis \iAijA2kB \Ib2) = K)ai|j)a2|A:)bi|0b2 gives 
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Pout 



^ Pijklmnst \iAljA2SBltB2){niAinA2kBllB2\ 
ijklrnnst 



Now, let us suppose vector 
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\4>rnn) = ~J= ( 1 0) Al \ m) A2 | 0) Bl \n) B2 

+|m) j41 |0) A2 |n) sl |0) B2 ) 
where m, n > and calculate the following mean value 



(0mn|Pout \<Pn 



-N 



Pm >r< 



(4) 



If Pm 7^ and r n =/= 0, this quantity is negative. On 
the other hand, the entanglement of discussed states can 
often be "masked" by the noise of original mixed states. 
E.g., conditional von Neumann entropy, S(p' A ) — S'(p), is 
positive for many particular cases here. Nevertheless, we 
will show that the entanglement is "strong" enough to 
violate CHSH-Bell inequality. 

The proposed conditional operation can be, in princi- 
ple, realized in the following way (see Fig. |l|) Let 
us assume a Mach-Zehnder (M-Z) interferometer with 
equal-length arms and with one photon in its input. Into 
both the arms of the interferometer we insert a nonlin- 
ear Kerr medium effectively described by the following 
interaction Hamiltonian 



H, 



Una) aa^ai 



(5) 



where and a are the creation and annihilation opera- 
tors of the mode corresponding to the left (or right) arm 
in the M-Z interferometer, a\ x and an, with i = A, B, 
are the creation and annihilation operators of modes Al 
(or Bl), and n is a real interaction constant. 
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FIG. 1. Preparation device; SA and SB denote the sources 
of thermal (pseudothermal) radiation, SPS is a single photon 
source, QND is quantum nondemolition measurement per- 
formed by the Kerr interaction and D+ and D_ are detectors. 

If there is a photon in the left arm of the central M- 
Z interferometer and the product KT; n t , where T mt is an 
effective interaction time,|j is set to be equal exactly to 



1 Quantity Ti nt has the meaning of the parameter of the de- 
vice. It is not a usual time variable. It represents the effective 
expression of the fact that the nonlinear medium has finite 
dimensions. 



7r then the described device realizes the phase shift ir in 
the left M-Z interferometer, A, and effectively flips the 
modes Al and A2 on the output. On the other hand, 
if there is no photon in the left arm then the states of 
modes Al and A2 stay unchanged, 



C/ j4 |n)Ai|0)A 2 |l) = |0)aiN.42|1), 
UA\n) A i\0) A 2\0) = |n)Ai|0)A 2 |0). 



(6) 



The same is true about the right arm of the central M- 
Z interferometer and modes Bl and B2. These unitary 
transformations Ui, i = A,B, can be expressed as 



Ui = U B s,iUldUBS,ii 



(7) 



where Ubs,i is the 50 : 50 beam splitter transformation 
and Ui t i accounts for the nonlinear interaction in Kerr 
medium, 



Ubs,i = exp 



n i rtrt ' 



Ui ti = exp(i7ra T aa^an) 



1^2) 



(8) 



So, if the photon goes through the left arm the modes 
Al and A2 are flipped while the state of system B is 
unchanged. Completely symmetrical situation occurs, if 
the photon goes through the right arm. 

Due to the path uncertainty of the photon in the in- 
terferometer the state of the whole system after the in- 
teraction is given by the formula 



|vf) = -^(|0>|1) \n) A i\Q)A2\0)Bi\m) B2 + 
+i|l)|0) |0)Ai|n) A2 |m) B1 |0)B2), 



(9) 



where the kets without any subscript denote possible 
states of the photon inside the M-Z interferometer sit- 
uated in the center. Which-way information is finally 
erased || by a beam splitter with amplitude reflectivity 
i/V2 (the last one in the M-Z interferometer) followed by 
two photodetectors D+ and D- (see Fig. [I]). Depending 
on which one of these two detectors fires we obtain one of 
two possible output states of modes Al, A2, Bl, and B2. 
Detector D + fires with probability w + = (1 + 5„o<5 m o)/2 
and if it clicks the following state is obtained 



|*+> = ^{\n)Al\0)A2\0)Bl\m) B 2 



-|0>Ai|n) A2 |m)m|0) 



B2- 



(10) 



Similarly, detector D- clicks with probability w_ = (1 — 
<5riO<5mo)/2 and when it fires one obtains the state 



!*_) = -^(|0)Ai|n) iU |m) B1 |0) 

-|n)Al|0)A 2 |0)Bl|m) B2 

which is exactly the considered state (0). 



B2, 



(11) 
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III. BELL-INEQUALITY VIOLATION 

In order to demonstrate the violation of Bell inequal- 
ities one needs local operations analogous to spin rota- 
tions. In our particular case the following operations do 
the job 

|n)i|0) a cos0 |n)i|0) a + sin0 \0)i\n) 2 for n + 0, 
|0>i|n) a ^ - sine |n)i|0) 2 + cosfl |0)i|n) a for n ^ 0, 
|0>i|0) 2 -|0)i|0) a , (12) 

where 6 is the parameter of transformation, it does not 
depend on n. 

Bell-type experiment consists of two "rotations" ac- 
cording to recipe JT^), performed by two possibly space- 
like separated observers, followed by realistic yes-no de- 
tection on each mode. Each such detection has only two 
possible outcomes (detector either fires or it does not), 
that can be described by projectors |0)(0| (for "no") and 

i - |o> <o | = Yln=i \ n )( n \ ( for "y es ")- Let us assi s n the 

following values to these outcomes: Zi=0 if the detector 
(in mode i) is quiet and Zj=l if it clicks. Then the results 
X and Y of local two-mode measurements (including "ro- 
tations" ) performed by the first and the second observer, 
respectively, can be expressed as 



= 0, 



X{0) = z AX {6) - z A2 (9), 
Y{9) = z m (9) - z B2 (9). 



(13) 



After the experiment is repeated many times and our two 
observers compare their results, the mean value of Bell 
operator (for CHSH inequalities) can be estimated, 

B = \C(6 A , B ) + C(6 A , 0' B ) + C{6' A ,6 B ) - C(6' A , # B )\, 

(14) 

where correlation function 

c{e 1 ,e 2 ) = Yj x * YkP{x 31 Y k \e A ,e B ) (15) 

3,k 

(summations go over all possible results). Every local- 
realistic theory || must fulfill the following inequality 
S < 2 ||. However, it follows from straightforward 
quantum-mechanical calculations that for state (^) the 
correlation function (|15f) reads 



c(e A ,e B ) = -cos[2(e A -e B )} (1 / o)(1 ro) . (ie) 

Therefore the results of the above mentioned local mea- 
surements performed on state (|^) violate inequality B < 2 
in principial. Maximal violation, 



B u 



2V2 



(l-p )(l-r ) 
1 - Poro 



(17) 
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(18) 



If both the mixed states have the same overlap with vac- 
uum state po = ro , the condition for the violation of Bell 
inequality for the considered angles is given in a simple 
form 



Po < 



V2-1 
V2 + 1 



0.1716. 



(19) 



As one can see the maximum value of B depends on the 
probability of the presence of the vacuum state in the 
input density matrices. Thus, if the input density ma- 
trices of systems Al and Bl do not contain the vacuum 
state the maximal violation of CHSH-Bell inequality is 
the same as for the pure EPR maximally entangled state 
of two spin-half particles. In the opposite case, the mean 
value of Bell operator decreases as the contribution of 
the vacuum state increases in the mixtures. It should 
be noticed that for properly chosen local measurements 
the violation of CHSH-Bell inequality does not depend 
on the randomness contained in the mixture but only on 
the overlaps of the vacuum state and the input density 
matrices. 



IV. THERMAL AND PSEUDO-THERMAL 
RADIATION 

There are two mixed states of special interest, namely 
thermal radiation, exhibiting Bose-Einstein statistics, 
and pseudo-thermal radiation, exhibiting Poissonian 
statistics. Let us study now the entangled states pre- 
pared by the device proposed in Sec. O when thermal 
and pseudo-thermal states are at the input. 

A single mode of thermal radiation has the density 
matrix 



in 



n 



^ (1 + (nW+r 

n v \ / / 



\n){n\, 



where 



(20) 



(21) 



For example, if the temperature of a radiation source 
(e.g., incandescent lamp) T ps 3000 K and the optical 
frequency u> « 2.5 x 10 15 Hz, the mean value of photon 
number is (n) ~ 1.77 x 10~ 3 . The probability of the 
vacuum state in the mixture is 



/ " ,= I " < ' M,( -^) = TT{ 



(22) 



occurs for the angles 



what leads to the value po w 0.9982 for the above given 
data. Thus in the optical region, the overlap of vacuum 
and thermal light is too large and the Bell-inequality vi- 
olation does not occur. 
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The dependence of the maximal Bell-inequality vio- 
lation on the parameter = ftaJi/ksTi, i = A, B of 
particular modes Al, Bl can be simply evaluated: 



2V2 



= 2V2 



exp(p A ) + exp(/3 B ) - 1 
1 



1 



(23) 



and it is displayed in Fig. g. Only for very small Pa 
and Pb, i-e., for high temperatures and small frequen- 
cies, CHSH-Bell inequality is violated. Thus for the 
given temperature T of both the thermal sources the in- 
frared component of radiation gives better results than 
the ultra-violate one. On the other hand, for the fixed 
frequency uj of both the sources the higher temperature 
leads to the stronger violation of Bell inequality. If both 
the sources are identical the Bell-inequality violation oc- 
curs only if the dimensionless parameter (3 satisfies re- 
lation (3 < In v/ ^ +1 ~ 0.1882 or the mean number (n) 

is sufficiently large (n) > 2(a/2 + 1) w 4.828. Conse- 
quently, for the visible component of radiation the ther- 
mal sources must have an "astronomical" temperature 
T > 101000K, whereas for the infrared component with 
w ~ 5 x 10 13 Hz, temperature T > 2021K is sufficient to 
obtain Bell-inequality violation. 




■<=> ^ 

FIG. 2. The maximal violation of the CHSH-Bell in- 
equality for thermal light as the function of parameters 
(3a = TiuJA/ksTA and (3b = TiuiB/ksTB 

Another interesting kind of mixed state is that corre- 
sponding to pseudo-thermal light fic| l . Its density matrix 
can be written as 



n) (n\ 



(24) 



mixture (|2J) is not vacuum state but it is state \n), where 
n corresponds approximately to the mean number of pho- 
tons (n) . Thus the overlap of pseudo-thermal light with 
the vacuum state is much less than for thermal light. The 
probability of the vacuum state in the density matrix ( [24] ) 
is po = cxp(— (n)). This leads to maximal Bell-inequality 
violation 



B = o /g [1 ~ exp(-(n) A )][l - exp(-(n) B )] 
1 - exp[-((n) A + (n) B )] 



(25) 



From Fig. || one can sec that in the case of pseudo- 
thermal light the Bell-inequality violation is achieved for 
less (n)A and (n)s than in the case of thermal light. If 
one considers two identical pseudothermal sources, then 
the Bell inequality is violated if (n) > In ^| + ^ . As the 
pseudo-thermal light with such a mean photon number 
can be experimentally achieved from the laser light in 
optical frequencies, the violation can be obtained more 
simply than for the thermal light. 



region of Bell-inequality violation 
for thermal light 



region of Bell-inequality violation 
for pseudo-thermal light 
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FIG. 3. The border of violation of CHSH-Bell inequality 
for thermal and pseudo-thermal light in dependence on mean 
photon numbers {n)A and (n)s. 

Real light sources emit to a large amount of different 
independent modes. The density matrix of this multi- 
mode state is given in the following form 



n 



00 



^ (1 + <^» 1+,V 



(26) 



where n M is photon number for particular mode [i and 
\iifj) is the Fock state of the corresponding mode. Let 
us suppose that this multi-mode thermal state is feeded 
to the inputs Al and Bl and the multi-mode vacuum 
states are present in the inputs A2 and B2^ The analysis 
presented in Sec. Ill may be generalized to multi-mode 



Pseudo-thermal radiation can be obtained from an 
intensity-stabilized single mode laser with the phase uni- 
formly distributed in the interval (0, 2ir). In contrast to 
thermal radiation, the maximally probable state in the 



2 In reality there could be a problem to set the proper pa- 
rameters of Kerr interaction for all the frequancy components 
together. 
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light in a straightforward way. We define the "rotations" 
of the multi-mode vacuum |{0}) and any excited multi- 
mode state |{n}) as follows, 

I W)i|{0}) 2 - cos9 |{n}) 1 |{0}) 2 + Binfl |{0}) 1 |{n}) 2 , 
|{0})i|{n}) 2 - - sine |{n}) 1 |{0}) 2 + cosfl |{0})i|{n}) a , 

(27) 

for {n} ^ 0, and for multi-mode vacuum in both the 
modes: |{0})i|{0}) 2 -> |{0})i|{0}) 2 . Detection that dis- 
criminates between the field vacuum and other states has 
two possible outcomes described by projectors |{0}) ({0}| 
and 1— 1{0}) ({0}|. It can be shown that the maximal vio- 
lation of Bell inequality exhibits the same form jl7| ) as in 
the case of single-mode radiation, but with the following 
notation 



Po 



n 



ro = 



n 



(28) 



With increasing number of the modes of thermal radi- 
ation the effective overlap of vacuum state and such a 
multi-mode field decreases and, consequently, the maxi- 
mal violation of Bell inequality is enhanced. In this way, 
the Bell inequality violation can be achieved for every 
thermal radiation, if the sufficient number of modes is 
taken into account. 



V. CONCLUSION 

An entangling device employing non-linear dynamics 
and postselection has been proposed and it has been 
shown that two mixed states can be entangled in such 
a way that the entanglement of the resulting state is 
strong enough to violate Bell inequalities (when proper 
local measurements are chosen). The disorder due to 
the statistical nature of the density matrices of input 
states is irrelevant - it does not influence the violation 
of Bell inequality. The only parameters affecting the 
maximum of the mean value of Bell operator are over- 
laps po = (0|p^i|0) and r = (0\pbi\0). This is also 
the reason of a contra-intuitive behavior when the en- 
tanglement increases as the input thermal state becomes 
more 'classical' ((3 — > 0), whereas in the 'quantum' limit 
(/? — > oo) the entanglement vanishes. Another contra- 
intuitive aspect of this phenomena appears if the multi- 
mode thermal radiation is considered. Since the overlap 
with multimode vacuum becomes smaller as the number 
of modes increases, the multi-mode thermal radiation can 
violate Bell inequality more notably, irrespective of its 
larger entropy. Thus this "classical-like" radiation can 
be strongly entangled in the ideal case and even exhibit 
the pronounced quantum nonlocality. Unfortunately, like 
the other kinds of mesoscopic states, the described quan- 
tum superpositions are very sensitive to the descructive 
influence of decohcrence and losses. 
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